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Abstract 

We discuss nonperturbative QCD evolution of nonsinglet nucleon structure functions, 
with particular application to the Gottfried sum. We show that the coupling of the quark 
partons to bound state mesons leads to nonperturbative contributions to the Altarelli- 
Parisi equations which, due to the axial anomaly, result in a strong scale dependence of 
nonsinglet structure functions for values of around the nucleon mass scale. We compute 
specifically the evolution of the first moment of the quark distribution, and find that it is 
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sufficient to explain recent experimental data which indicate a violation of the Gottfried 
sum rule. 
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1. Violation of the Gottfried Sum Rule 

In the wake of the furore over the EMC measurement |^ of quark polarized structure 
functions, experimental results from the NMC0 for the Gottfried sum again appear to be 
in contradiction with theoretical expectations. More specifically, the NMC measure the 
ratio of the cross-sections for unpolarized deep-inelastic scattering from deuterium and 
hydrogen targets. The ratio of structure functions -F2 (x, Q^) / -Ff (x, Q^) is extracted in the 
range 0.004 <a;<0.8at(5^ = 4 GeV^. Extrapolating to small x using Regge behaviour, 
and assuming further that the neutron structure function is just F2 = — F2 allows a 
determination of the Gottfried sum, defined as 




with the result 

5g = 0.24 ±0.016 at Q^=4GeY^. (1.2) 

This result (unlike previous less precise experimental determinations which had a consistent 
central value but larger errors) is in striking contradiction with the Gottfried sum rule, 
which would have Sq = ^ 0- 

Both of the assumptions used by the NMC to arrive at the result ( |1.2| ) may be ques- 
tioned: on the one hand, if Regge behaviour is assumed to set in only at very small x 
(~ 10""^ or less) it is not difficult to increase Sq [Q, while on the other if nuclear shad- 
owing is taken into account the neutron structure function is certainly larger than the 
difference between deuterium and proton structure functions, thus leading to a yet smaller 
value of Sg, perhaps^] around or even below Sq ~ 0.2. The settling of these issues rests 
ultimately with experiment, and in particular with the availability of better data in the 
small-x region. However, if we accept that the result ( |1 . 2| ) is correct, perhaps with slightly 
underestimated error, we must conclude that the Gottfried sum rule does not hold at the 
scale of the NMC measurement. 

Let us therefore review how this sum rule is derived. In the parton model, F2{x) 
is related to the sum of parton densities, weighted by the electric charges squared (see 
e.g.§): 

F2ix) = xj^ej {q,ix) + q,ix)) , (1.3) 

i 

where qi{x) (qi(x)) is the quark (antiquark) density of flavor i. This identification is exact 
(by definition) to all orders in the QCD parton model. Although the first moments of qi{x) 
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and qi{x) diverge because of their small-x behaviour, this divergence is expected to cancel 
in nonsinglet combinations of quark distributions (such as in the Gottfried sum), on the 
basis the smaU-x behaviour of qi{x) expected from Regge theory 0. 

It foUows that the Gottfried sum measures the difference between the sums of the 
square charges of partons (quarks plus antiquarks) in proton and neutron. Assuming exact 
isospin symmetry 

Sg = dx {el-ea){Qu{x) -Qdix) +qu{x) -Qdix)) = l{qi + qi), (1.4) 

^0 

where qi and qi are the first moments of the nonsinglet quark distributions 

q{x) = quix) - qdix), q{x) = qu{x) - qd{x) . (1.5) 

On the other hand, one may derive an exact sum rule (the Adler sum rule) for the charge- 
conjugation odd combination of first moments which counts quarks minus antiquarks. By 
charge conservation this is fixed: 

SA = qi-~qi = 1- (1-6) 

Only valence quarks contribute to the Adler sum rule, since (by definition) q^^ = 
while the Gottfried sum receives contributions from both valence and sea quarks. However 
if we assume further that the quark-antiquark sea is isotopically neutral, then = qf^^ = 
0, and we find (using ( pTBI) in (|1.4| )) the naive result Sq = ^ alluded to above. 

The NMC result (|1.2| ) thus suggests that the light fiavor content of the nucleon is very 
different from the standard folklore Q of an asymmetric valence quark content, plus sea 
quark, antiquark and gluon distributions which respect all possible symmetries. In par- 
ticular, a detailed analysis]^] indicates that isospin violations (such that the total number 
of sea quark-antiquark pairs in the proton and in the neutron are not exactly the same) 
cannot by themselves account for the discrepancy; the value ( |1.2| ) requires a large violation 
of the SU(2) symmetry of the quark-antiquark sea, such that the number of up quark- 
antiquark pairs in each nucleon must be different from the number of down pairs. More 
precisely, the NMC data require ||] the isospin violation to be ^ 5%, and the sea asym- 
metry to be opposite in sign to the valence asymmetry, and about 30% of its magnitude. 
Thus (recalling eqn.(|1.6|)) qf^^ + qf^^ ~ —0.3, and not zero as naively assumed. 
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This conclusion is not only unsettling from the phenomenological point of view, given 
that the naive assumption of a symmetric sea is embodied in all phenomenological param- 
eterizations of parton distributions which are used in analyzing deep-inelastic scattering 
dataS, but is also puzzling from the theoretical viewpoint. Indeed, one might expect the 
quark-antiquark sea probed in deep-inelastic scattering to be mostly produced through 
QCD evolution. But the scale dependence of Sg as computed perturbatively at two loops 
in QCD[^ is negligibly weaki, leading to a variation of Sq by less than 1% in the per- 



turbative region. Furthermore, the same computation ||1 1|| shows explicitly that the naive 
explanation ||12|| of the sea asymmetry as a consequence of the Pauli principle (according 
to which the emission of sea up pairs in the proton should be reduced due to the up excess 
in the valence component [jl^) is incorrect: due to subtle interference effects the Pauli 
antisymmetrization actually leads to a slight perturbative enhancement of the sea flavor 
asymmetry with the same sign as the valence asymmetry. 

It has been suggested [|T2|,|T^,^ that the effect may instead be related to the presence 
of a pion cloud in the nucleon wave function. Indeed, if one assumes that some portion 
of the sea is produced through transitions where a pion is radiated by a nucleon, then a 
proton would seem to favor the transition where a tt~^ is created, namely p ^ n + tt"*", 
over that where a tc~ is created, namely p A~^~^ + n~ (and conversely for a neutron), 
because of the nucleon-A mass difference. This, spelling out the quark content, is seen to 
favor the production of dd pairs over uu pairs, thus producing a flavor asymmetry with 
the right sign. The observed effect may also be accommodated rather easily in most of 



the various effective models of the nucleon; examples are the Skyrme model||^,16|, chiral 
soliton models and various bag modelspl JT^] . 

The problem with this kind of explanation is not only that it is difficult to make 
quantitative and model-independent predictions, but that the connection between the 
models (which are presumably only valid at some unknown but low energy scale) and the 
parton distributions measured by the NMC is not at all clear. Furthermore, even accepting 
that pion effects are somehow responsible for the observed asymmetry, it is not enough 
to assume that the tt"'" and 7r~ clouds in a nucleon are unequal, since each pion gives a 
vanishing contribution to the flavor asymmetry of the quark sea. Rather, one has to make 
a separate dynamical assumption on the pion production mechanism, in order to obtain 



^ Recently new parameterizations have appeared, which incorporate the NMC asymmetry [ 10 1. 
^ Though not zero, as seems to be sometimes believed (see for example ref. Q). 
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an anticorrelation between valence and sea asymmetries. Such a mechanism has been 
presented in ref. []T^, where it is shown that if (part of) the sea is produced by dissociation 
of a quark into a quark plus a pion, then the breaking of the axial U(l) symmetry in the 
pseudoscalar spectrum favors flavor-changing transitions over flavor-preserving ones. A 



calculation [|T9[ based on the chiral quark model of ref.[^ seems to support this idea. 

It is interesting to compare this situation with a similar instance of the failure of 
naive expectations, namely that of the EUis-Jaffe sum as measured by the EMC|]l|] (for 
reviews see ref.|^T|). In that case too the first moment of a nucleon structure function 
(the polarized proton structure function gf) is found to be considerably smaller than the 
parton model expectation, implying that the quark contribution to the "proton spin" (more 
precisely, its isosinglet axial charge) is consistent with zero. There, however, perturbative 
QCD suggests a partial explanation of the observed effect: due to the axial anomaly, the 
isosinglet axial charge evolves, despite being classically conserved. This evolution may 
be interpreted in the parton model by writing the charge as a sum of a polarized quark 
contribution (which does not evolve) and an anomalous gluon contribution (which does). 
Whereas the perturbative evolution (which, again, only occurs at two loops) is too weak to 
explain why the gluon contribution is so large, the observed discrepancy could be explained 
by a nonperturbative dynamical mechanism which generates a gluon contribution of the 
appropriate sign and magnitude. Such a mechanism would incorporate the same infrared 
non-perturbative physics that resolves the U(l) problem by giving a mass to the flavor 
singlet pseudoscalar meson (the rjo) [^. It could take the form of infrared vacuum effects 
(as modelled for example by instantons) which reduce the starting value of the charge. 



or a non-perturbatively induced scale dependence of the isosinglet axial charge |2J] . 

Here we will show that the anomaly also affects the non-perturbative evolution of 
the Gottfried sum, albeit in a rather more subtle way. We compute the non-perturbative 
scale dependence of nonsinglet quark structure functions by including bound state emis- 
sion, thereby generalizing the usual Altarelli-Parisi evolution equations. We find that 
when the axial anomaly is taken into account through the generation of a mass for the r]o 
(as suggested in ref. fl^), then a flavor asymmetry in the quark sea is generated dynami- 
cally and the Gottfried sum acquires a rather strong scale dependence over intermediate 
scales. At very large and very small scales this dependence flattens out, so that the two 
loop perturbative evolution eventually takes over at large Q^, whereas a constant "quark 
model" value is attained in the infrared region. We find that if this is identified with the 
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naive expectation Sq = 1/3, then the experimental value of Sq, eqn .( 11.21) , can be easily 
reproduced. 

The paper is organized as follows: in section 2 we review the general theory of the 
dependence of the nonsinglet structure function F2, both from the point of view of the 
operator-product expansion, and the QCD parton model; in section 3 we generalize the 
latter approach to include non-perturbative coupling to mesonic bound states, we describe 
the determination of the nonperturbative contribution to the anomalous dimensions, and 
we discuss how this can lead to flavor symmetry breaking effects due to the U(l) anomaly; 
in section 4 we compute the non-perturbative contribution to the appropriate splitting 
functions by explicit determination of the relevant cross-sections; in section 5 we determine 
the dependence of the first moment of the quark distribution, hence of the Gottfried 
sum, and compare our results with the NMC data; and finally in section 6 we summarize 
our results and discuss how this approach may be developed into a more general technique 
for the computation of nonperturbative contributions to structure function evolution. 



2. Evolution of Nonsinglet Quark Distributions 

Because the structure functions F2{x) are proportional ( |1.3| ) to the quark distributions 
qi{x), the problem of the determination of their dependence reduces to the computation 
of the QCD evolution of the latter. Since we are only interested in the difference — 
we need only discuss the evolution of the nonsinglet distributions q{x) defined in 
( |T75|) , although most of what follows would apply equally to other nonsinglet distributions. 

We begin with a discussion of the evolution as dictated by the operator product ex- 



pansion |jTT|]. This is not straightforward, because the (electromagnetic) structure functions 
F2 are even under charge conjugation, as eqn.( p..3| ) demonstrates explicitly. It follows that 
at leading twist only the even moments of F2 can be expressed as the matrix element of a 
local operator, since the leading-twist contribution to the A^-th moment has spin hence 
charge conjugation C = {—)^ . As shown in ref.|ll], however, this does not imply that it 



is impossible to give a meaning to odd moments of F2 or their evolution, as is sometimes 
claimed [|2^ . 

Indeed, consider the combinations 

q^{x;t) = q{x;t)±q{x;t), (2.1) 
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where we have indicated exphcitly the dependence on the scale t = InQ"^ j (j?. Since are 
eigenstates of charge conjugation, they evolve independently. Even moments of and 
odd moments q~ have the scale dependence 

q%{t) = AAr(t, to) g^(io), N even, 

(2.2) 

QNi^) = ^N{i,to)q]si{to), N odd, 

where q^{t) = dx x^~^q^{x;t) is the A^-th moment of q'^(x;t), to = InQo/Ai^, and the 
scaling factor A]\[{t, to) is given in terms of the anomalous dimensions 7Ar(t) of the twist-2, 
spin-A^ nonsinglet local operators by 

lnA^(t,to) = / 7N{t')dt'. (2.3) 

Jto 

If we let 

7A^ = 7^^ + (-)^7^^ (2.4) 

and assume that both 7^ and 7^ admit an analytic continuation for all values of (as may 
be verified explicitly in perturbation theory [|TT| ), it follows that the evolution equations 
( p.2|) also admit a continuation for all values of A^, such that q^ evolve multiplicatively 
according to 

4(t) = A±(t,to)4(^o), (2.5) 
where A^(t,to) are defined in terms of the anomalous dimensions 

7^W=7^^W±7^'(t)- (2.6) 
In particiilar, the Gottfried sum evolves multiplicatively as 

SG{Q^)^At{t,to)SG{Ql). (2.7) 

This somewhat formal result becomes more transparent in the parton model, where 
the evolution of the quark distributions is found by using the Altarelli-Parisi equations 

j j 
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where V are splitting functions, the sum runs over flavors, and V ® q = J^^'P{^)<liy)- 
Since, by charge conjugation, Vq^q. = Vq^q■ and Vq-q. = Vq^q., if we assume all relevant 
flavors to be massless there are only four distinct quark-quark or quark-antiquark splitting 
functions, namely V^q,Vqq,Vqq^ and Vqq^ , where (p^^^ is any Vij such that i = j 
{i 7^ j). The evolution of the nonsinglet quark and antiquark distributions ( p. .51 ) is then 
simply given by 



d 

Jt 



— q = Qqq 'S)q+ Qqq ® <?, 



^q=Qqq®q+ Qqq ® Q', ^^"^^ 

This shows immediately that the eigenstates of the Altarelli-Parisi evolution are g^. 



J^q^ = iQqq±Qqq)^q^. (2.10) 



eqn. (|2.1D , since 



Taking moments of this equation and comparing the result with ( |2.5| )-( PTB| ) it is clear that 
the anomalous dimensions 7^ ( |2.4| ) are just the A^-th moments of the evolution kernel on 
the right hand side of eqn.( p.lOD , whereas 7^ (In) N-th moments of Qqq (Qqq); 

this may be confirmed in perturbation theory by explicit computation. In particular, 
the nonsinglet evolution of F2 (and hence that of the Gottfried sum) is found by taking 
moments of Qqq + Qqq. 

In the parton model several useful perturbative properties of the anomalous dimen- 
sions 7^ become immediately apparent. Firstly, at one loop there is only one diagram 
which may contribute to nonsinglet evolution, namely gluon radiation by a quark line. 
Hence all the splitting functions Vqq vanish. On the other hand, at two loops the split- 
ting functions and do not vanish, since a quark may radiate a gluon which in 
turns radiates a qq pair (see fig. 1). Furthermore the splitting functions V^q and V^q^ 
are in general different already at two loops: the cross sections for fiavor-diagonal and 
fiavor-non-diagonal production differ since in the flavor diagonal case there are two iden- 
tical particles in the flnal state, which must be antisymmetrized. It follows that all the 
7^ start at two loops. Furthermore, specializing to the flrst moment, charge conservation 
implies that the flrst moment of q~ must be conserved (and thus that the Adler sum ( |1.6| ) 
does not evolve). It follows that 7^"^ = 7^'^, and in particular that they both vanish at one 
loop. 
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In conclusion, we see that in perturbative QCD the evolution of the Gottfried sum 
is governed by the anomalous dimension 7^ = 27^"^, which is the first moment of the 
splitting functions for diagonal minus non-diagonal emission. This quantity starts at two 
loops, thus leading to a perturbative contribution to the evolution factor of the form 

(A+)pert(t, to) = l + -^{asit) - as{to)) +■■■ (2.11) 

where 7^*^^'* is the two-loop coefficient in 7^*", and b^^^ is the one-loop coefficient in the 
beta function. Explicit computation shows that ||l 1|| "f^^'^yb^^^ ~ 0.01; as mentioned in the 
introduction the perturbative evolution is thus entirely negligible unless as is so large that 
perturbation theory is useless. 

This result was derived in the flavor-symmetric, vanishing quark mass limit. In prin- 
ciple, it may be corrected by two different kinds of effects. Firstly, if current quark masses 
are taken into account, then the simple evolution equation ( |2.5D is corrected by flavor- 
violating terms. These are usually assumed to be negligible, in that mass effects should be 
suppressed by powers of m^/Q^, although they have never been determined explicitly. If 
one also admits the possibility of isospin violation, then Sq also contains a term propor- 
tional to the proton-neutron difference of singlet structure functions, and the evolution 
of these must also be taken into account; Sc will then no longer evolve multiplicatively. 
This effect has been discussed in ref.[0, and is negligible for realistic amounts of isospin 
violation. We will neglect both these effects henceforth and assume isospin symmetry to 
be exact. 



3. Evolution by Bound State Emission 

The physical reason for the smallness of the perturbative evolution found in the previ- 
ous section should be readily apparent: if we were to neglect final-state antisymmetrization 
effects, then the interaction which governs the QCD evolution would be fully symmetric un- 
der the U{Nf) fiavor group, and the Wigner-Eckart theorem would imply that = . 
There would then be no nonsinglet evolution, and the dynamically generated sea would 
be fiavor symmetric, as naively expected. Final-state antisymmetrization effects spoil this 
argument, and provide a (tiny) dynamical violation of the flavor symmetry. 

However, a much larger violation appears if we take into account ||19|| that in QCD the 
U{Nf) flavor symmetry is dynamically broken due to the axial anomaly ||22|| . In particular, 
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consider diagrams where a quark couples directly to a qq bound state (fig. 2). We may view 
these as generated dynamically from diagrams which have the structure of the perturbative 
ones of fig. 1, but where the emitted (and unobserved) quark-antiquark pair may interact 
non-perturbatively. These diagrams will be suppressed at large scales where, due to the 
extended nature of the emitted bound states, they must lead to higher-twist contributions 
to the evolution equations — contributions suppressed by powers of l/Q^ in the large 
limit. However at low enough scales their contribution to nonsinglet evolution will 
probably be much greater than the perturbative diagrams. At very low scales it is no 
longer clear that the Altarelli-Parisi evolution picture is applicable at all. There might 
nevertheless be an intermediate region where the Altarelli-Parisi picture applies and the 
evolution due to bound-state emission is significant. 

First, consider the situation in the chiral limit, with Nf massless quarks. Then, it is 
clear that all of the bound state emission diagrams will respect U{Nf) fiavor symmetry, 
and thus will not contribute to nonsinglet evolution, with the sole exception of those in 
which pseudoscalar mesons are emitted. This is because it is only in the spectrum of 
pseudoscalar mesons that U{Nf) is broken down to SU{Nf) by non-perturbative effects 
coupled through the axial anomaly [p2| ; the Nj — 1 Goldstone bosons will make a significant 
contribution to nonsinglet evolution which is not cancelled off by the contribution of the 
singlet pseudoscalar (the ryo) since this is kinematically suppressed. Away from the chiral 
limit, the emission of other mesons will also begin to contribute, but their contributions will 
remain relatively small mainly because the mass splittings are themselves relatively small, 
but also because they are in any case kinematically suppressed when compared to those of 
the light pseudoscalars. It follows that if we concentrate on nonsinglet evolution, then of 
the infinite tower of diagrams for bound-state emission (fig. 2) only those corresponding 
to pseudoscalar meson emission need be considered. 

To make this discussion more quantitative consider the inclusion of pseudoscalar bound 
states in the evolution equation ( |2.8| ). We now have (omitting gluon contributions for 
simplicity) 



d 




3 



a 




(3.1) 



d 
dt 



a 
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where 11" (x; t) are the distribution functions of the pseudoscalar mesons (the pion octet and 
the r]'). The first two equations in ( |3.1| ) express the evolution of the quark and antiquark 
distribution due to pseudoscalar emission, whereas the last equation gives the evolution of 
the pseudoscalar meson distribution. The former is determined by the splitting functions 
Vqq, which expresses the probability of a quark to be emitted by another quark (with 
emission of an unobserved pseudoscalar), and "P^n, which expresses the probability of a 
pseudoscalar to fragment into a quark and an antiquark (one of which is observed). The 
latter is determined by the splitting functions Vuq and Pnn, which give the probability of 
a pseudoscalar to be emitted by a quark or another pseudoscalar, respectively. 

Isospin and charge conjugation invariance imply that for neutral pseudoscalars 11^ = 
(7r°, ?7, V) 

'Puno = 'Pdiio = 'PuW = '^dno; (3-2) 

whereas for charged mesons 

'Pu7r+ T^diT^ "^^77+ T^u-K" — Pqin 

(3.3) 

^d7r+ Pun~ PuTT+ 'Pd-K~ ^" 



The non-perturbative evolution of the nonsinglet quark and antiquark distributions ( pTq ) 
and the pion distribution 

7r(a;; t) = 7r''"(a;; t) — 7r~(a;; t), (3-4) 



is thus given by 



d 

It 
d 

lit 
d 

~dt 



— q=Qqq'S)q+Qqq(S)q+ Vq^ ® TT, 

^q = Qqq(S)q + Qqq(S)q~Vq^(^n, (3.5) 

TT = Pvrg ® (? " ^) + Ptttt ® TT. 



The combinations ( |2.1| ) then satisfy 

d 



= {Qqq+Qqq)®q'^, 
J^q~ = {Qqq- Qqq)®q~ +Vq^®2Tl, (3.6) 

d ^ - ^ 

^71" = A'tt? ®q + A'tttt ® TT, 
11 



where, to first order, only the sphtting functions Vqq = Vqq, Vqu, Vqu, Vuq and Vuq are 
nonvanishingi (and in particular Qqq = 0). 

To actually compute the nonsinglet evolution due to pseudoscalar meson emission, 
we use the probabilistic interpretation[^ of the Altarelli-Parisi equations in which the 
splitting functions Vqq are obtained through the same sort of procedure which leads to the 
perturbative splitting functions]^. That is, QCD evolution is in general due to the fact 
that in deep inelastic scattering a parton may radiate another particle (perturbatively a 
gluon, but for present purposes a pseudoscalar bound state, H) either before or after in is 
struck by the virtual photon. The splitting function Vq,.q.{x) expresses the probability of 
the quark parton qj to carry the fraction x momentum of the quark parton qi after the 
radiation process, and it may be obtained obtained from the cross-section for Compton- 
like scattering of a (virtual) photon 7* before or after emission of a bound state (fig. 3). 
This corresponds to an effective resummation of the same ladder diagrams which lead to 
the usual Altarelli-Parisi evolution, with the exchanged gluons replaced by pseudoscalar 
mesons. These diagrams may be shown to give the dominant contribution in the leading 
logarithmic approximation in a theory where quarks are coupled to pointlike pseudoscalar 
particles in just the same way as for perturbative QCD, but without the complications 
arising from gauge invariance [0,^ . In section 4 we will discuss the applicability of this 
picture to extended pseudoscalars. 

The contribution of the emission of a bound state 11 to the splitting function Vq^q. is 
then given by 

[^....(^;^)]n = |<^"(^;^)' (3-7) 

where a'^.^{x;t) is the total cross section for emission of the state 11 with quantum num- 
bers corresponding to the given splitting function, expressed in terms of the usual scaling 
variables and integrated over all k±. This cross-section is adimensional, due to the ex- 
traction of a scale factor 0. Notice that the so-called "loss" terms [|6|,p6[|, namely those 
terms which give the decrease in the probability of finding a quark with momentum frac- 
tion X due to the fact that if such a quark undergoes a radiation process it is lost from 
the observed momentum interval, are absent from eqn.( |3.7| ), which gives the full splitting 



^ "Pnn begins at second order (just as the usual perturbative contributions to Vqq only begin 
at second order), while although Vqq may receive nontrivial contributions from diquark emission, 
the U {Nf) symmetry of the diquark is not broken anomalously, and thus they may be ignored for 
the same reason as we ignore the contribution of all mesons other than pseudoscalars to Vqq. 
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function. This is because when a meson is emitted the quark which is struck by the virtual 
photon is never the same as the original one, even if it has the same flavor. 

The non-perturbative contribution to the anomalous dimension 7^'^'' is then found 
by taking moments of the splitting function: 

/* ^ m ax ^ P^m ax 

bTit)]n = dxx"^-' [P,,,,(a:;i)]n = jj^ dxx^'-'al^ix-t). (3.8) 

Notice that the upper limit of the ^-integration is not 1 but rather 



1+ n2 



M2 



where M is the mass of the pseudoscalar bound state. The cross-section vanishes for 
X > a^max- By continuity we also expect cr'^.^ix^ax'-i't) = 0; this justifies the exchange in 
the order of differentiation and integration in the second expression ( |3.8D . The anoma- 
lous dimensions 7^ and 7^"^ are similarly determined by taking moments of the splitting 
functions Vqn and Vjiq obtained by differentiating the cross-sections o"gn and anq- 

Besides the usual crossing symmetry relations among splitting functions 0, charge 
conservation as expressed by the Adler sum rule eqn.(L6) imposes a nontrivial consistency 



condition on the first moments of the various splitting functions. Combining with 
the evolution equation ( |3.6|) and taking the first moment gives 

7r(t) + 27rW^iW = 0, (3.10) 

where 7^^*^ is the first moment of Qqq, and qi{t) and 7ri(t) are the first moments of q{x;t) 
and Tc{x; t) respectively. If we now assume that there exists a scale to such that ni{to) = 0, 
the entire pion distribution 7r(x; t) is generated dynamically through the relevant evolution 
equation ( |3.6|) . Integrating up the first moment of this evolution, and using the result in 
eqn.( |3.lOD leads to the condition 



7r(t) + 27r(t) H,{t) - <(to)] = 0, (3.11) 

where o"^^ is the first moment of the appropriate cross-section, so that 7f'^(t) = ^c^g- 
Now eqn. (|3.11| ) means that 7i^(io) = 0; the assumption that vri(t) is entirely generated by 



the evolution eqn. (|3.6D is consistent only if the evolution of qi{t) flattens at small t. The 



scale to at which the evolution flattens is then the same as that where 7ri(t) vanishes. 
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In this formalism the so-caUed SuUivan process[^|, where the virtual photon scatters 
inelastically on the quarks in a pion which is in turn carrying a fraction of the nucleon 
momentum, is regarded as a measurement of the n"(a;) structure functions through their 
fragmentation into quark-antiquark pairs, as described by the splitting function Vqu- This 
is the nonperturbative analogue, for pseudoscalar bound states, of a perturbative measure- 
ment of the gluon distribution via photon-gluon fusion. Because some of the bound states 
we consider are charged, there is also the possibility of a direct coupling of the virtual 
photon to these states, corresponding to processes where the pion does not fragment inco- 
herently, but rather is viewed as a pointlike effective constituent. Such processes, however, 
are suppressed by the square of the pion form factor, evaluated at the photon's virtuality 
Q^, and are thus negligible |]2^. Furthermore they are only relevant to the evolution of 
singlet structure functions; for nonsinglets, and in particular for the Gottfried sum, they 
vanish identically because of isospin. 

Eqns (|3.1|) - (|3.7| ) provide us thus with a well-defined framework in which to calculate 
nonperturbative contributions to the evolution of structure functions due to the inclu- 
sion of pseudoscalar mesons as effective constituents in the ladder diagrams which lead 
to Altarelli-Parisi evolution. In order to avoid double counting, the non-perturbative 
contribution to the anomalous dimensions determined in this way should not be added 
directly to the usual perturbative anomalous dimensions. Rather, the non-perturbative 
contribution, while dominant at small Q^, must fall as a power at large Q^, eventually 
falling below the logarithmic perturbative evolution. At around this point (which would 
ideally be in a region where both calculations can still be believed) the non-perturbative 
and perturbative anomalous dimensions should be matched together. In the particular 
case of the Gottfried sum the perturbative evolution eqn.( |2.1I| ) is so tiny that it may be 
ignored; the only significant evolution is nonperturbative. 

We specialize now to the evolution of q^{x;t). It is apparent from eqn.( |3.6|) that 
the pion distribution n{x; t) makes no direct contribution to this evolution. This is to 
be expected on the basis of a simple quark counting argument: because quarks and 
antiquarks contribute with the same sign to Gottfried sum eqn. (|l.lj ), (|1.3|) , it follows that 



any pion always give a vanishing contribution to Sq eqn.(Ll). The emission of bound states 
contributes nevertheless indirectly, by leading to nonperturbative evolution of {x;t). The 
evolution of the Gottfried sum discussed in the present approach thus appears to be due to 
a mechanism which is physically distinct from the Sullivan process; indeed, in the present 
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approach scattering on the pion component of the nucleon does not contribute at aU to 



Sg, contrary to the suggestion in ref.[|15 



The non-perturbative evolution of the A^-th moment of the nonsinglet quark distri- 
butions q^{x) ( ]1.5D ,(P7TD is found by computing the anomalous dimensions ( ^.81 ) for the u 
and d flavors and the pseudoscalar mesons tt^, tt^, t) and ry'. For simplicity we may assume 
that the latter are pure octet and pure singlet, since the appropriate mixing angle 6'p is 
small (sin^^p ~ 0.1). Noting that the probability for flavor-diagonal emission by a light 
quark is obtained by adding the probabilities for the process to go through rr^, ry and ry' 
meson emission as dictated by their quark content we find 
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where aj^ ^ is the A^-th moment (defined as in the right hand side of eqn.( |3.8| )) of the cross 
section for the process of fig. 3, and we have included for completeness strange quarks and 
kaon emission. Using eqn.(^7T^) in the evolution equation ( |3.6D we finally get 

Tat — Tat — In In — 

d 

^ Jt U-A^ 

The non-perturbative contribution to the evolution factor for the Gottfried sum (|2.7|) is 
thus 
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(3.14) 



^ Notice however that in some specific low energy models the Sullivan process is actually 
related|14] by the dynamics of the pion radiation mechanism to the distribution of quarks in the 



recoil hadron after pion emission, thus it may contribute indirectly to quark distributions in a 
way which is analogous to that which we discuss here. In these models, however, there is no scale 
dependence of structure functions and quark distributions. 
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Whereas the perturbative evolution ( p.ll ) was very smaU the non-perturbative scale 
factor (|3.14| ) is potentially large. In particular, the role played by the anomalous U{Nf) 
symmetry breaking in the pseudoscalar spectrum is now apparent: because of the large 
ry-TT and r]'-n mass differences we expect the corresponding emission cross sections to be 
significantly different. In the t/(A^/) -symmetric limit all cross sections would be the same 
and the right hand side of eqn. (|3.13D , and thus the exponent in eqn.( p.l^ , would be 
identically equal to zero. In reality, whereas for very large values of the effects of the 
pseudoscalar mesons' masses should be negligible and the non-perturbative scale factor 
should flatten (as l/Q^), for scales of the order of the nucleon mass the effects of the mass 
splitting will be so large that the r] and i]' production cross sections will be negligible when 
compared to that of the tt. Thus because of the axial anomaly the diagonal radiation 
process is disfavoured compared to the non-diagonal one, as half of it proceeds through 
r] and rj' emission, which is dynamically suppressed. So not only is the evolution purely 
multiplicative, but since the cross-sections are all positive, it necessarily results in the 
screening of the Gottfried sum; eqn.( |3.l4D shows that Sq is always reduced as is 
increased. 

In order to make these observations quantitative, we must proceed to the computation 
of the various radiation cross sections and the determination of their flrst moment. We 
will do this in the next section. 



4. Non-perturbative Splitting Functions 

The computation of the non-perturbative splitting function reduces to the compu- 
tation of the cross section for the process 'y*q — > Uq, where 11 is a pseudoscalar meson, 
which at leading order is given by the two diagrams displayed in flg. 3. In order to perform 
this computation we must introduce an effective quark-pseudoscalar coupling. Writing the 
quark propagator as S~^{p) = + E(p^))~^, where is the quark self energy, the 

most general such coupling is, by definition, 

X{k,p) ^ S-\p + lk){Omp - lk)i;{p + lk)\U{k))S-\p - Ik), (4.1) 

where |n(/c)) is the meson state, and the momenta are notated as in fig. 2. The quark 
dynamics may thus be summarized by writing an effective action 

5eff = jd''pi;{p)i]pij{p) + jd'^p jd''ki;{p+^k)fx{k,p)U{k)i;{p-lk), (4.2) 
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where U{k) = exp (i'j^Il{k) / fu) is a nonlinear representation of the pseudoscalar field 
n(/c); expanding U in powers of 11 gives the quark self energy and then couplings of the 
qq pair to increasing numbers of pseudoscalars H. In the chiral limit (i.e. ignoring quark 
masses) the quark self-energy and the on-shell coupling to pions are thus related by the 
chiral Ward-Takahashi identity 

Sb) = /x(0,p), (4.3) 



/ being the meson decay constant (see for example [^^). Effective actions such as eqn 



have long been used to discuss dynamical chiral symmetry breaking and to compute meson 
form factors and decay amplitudes [|3T|-|33| ; more recently attempts have been made to find 
a formal justification for them within QCD ([^ and references therein) and to develop 



their low energy phenomenological implications more systematically ||35|| . 

The effective coupling (Of) may be expanded in terms of a set of four scalar vertex 



functions (f)i{k,p), i = 0, 1, 2, 3; 

X{k,p) = Mk,p) + Mk.pW + Mk,p){k ■ p)i> + (/)3(/c,p)i(# - ^|^). (4.4) 

The vertex functions are all even functions oi k ■ p, and for on-shell pseudoscalars we may 
expand them as 

oo 

(/>.(fc,p) = 5^</>^(p2)(fc■p)2^ (4.5) 

n=0 

Now, in principle, all of the on-shell vertex functions (p^{p) may be computed by solv- 
ing the homogeneous bound state equation, separated out into a hierarchy of equations 
corresponding to increasing powers of k. In practice of course this is not possible since 
the quark-antiquark scattering kernel is unknown except at large Euclidean momentum 
transfer. The asymptotic behaviour for large Euclidean p^ is given by ||36|| 



(logp^)'^o 
pCoo (p2)n+l ' 



(4.6) 



— * oc 



where the numbers Sf are related to the anomalous dimensions of various local operators, 
and decrease monotonically in n. Asymptotically, the dominant vertices are thus (pi{p'^) = 
4>^{p'^) for i = 0,1,2. Furthermore, it can be seen from examination of the bound state 
equations that all the <^'^{p^) remain finite as ^ 0. 
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The normalization of the vertex functions is fixed by the chiral Ward-Takahashi iden- 
tity (^4.3[ ) in terms of the self-energy and /^r: 



^oip') = np')/U- (4.7) 

The scale of S(p^), in turn, is in principle given by solution of the (nonlinear) Schwinger- 
Dyson equation in terms of Aqcd- Infrared uncertainties make such a determination 
unrealistic, however [p6|, and it is more practical to use the Pagels-Stokar condition[p3| 
which normalizes S(p^) to the pion decay constant /^r = 93 MeV: 

12 rfS(p2) 



ilr^ Jo ^ ^ ' (p2+E(p2)2)2 



Then, using eqn. ( [4.7|) in eqn. (|4.8| ) fixes the strength of the meson-quark coupling, or (more 



precisely) the normalization of the vertex function in terms of the scale of the quark self- 



energy. A consistency check on the condition (|4.8[ ) can be obtained by considering the 



asymptotic behaviour of the pion electromagnetic form factor (also derived in ref . ) : 

In principle all the other vertex functions may now be found using the homogeneous Bethe- 
Salpeter equations without the need for further normalization conditions, though in prac- 
tice this is again not really viable since they are too sensitive to the unknown infrared 
dynamics. 

It is now straightforward to compute the spin-averaged cross section corresponding 
to the two diagrams of fig. 3. Assuming that all the quark partons are massless, writing 
M for the mass of the emitted meson, and —Q^ for the virtuality of the incident photon, 
the usual Mandelstam invariants s, t, u, are given, in terms of the four external particles' 
four-momenta indicated in fig. 3, by 

s = {p^ + qf = {Pf + kf = -Q2 + ■ q = M2 2pf ■ k, 

i = {p, - kf = {pf - qf = -2pfq = M^- 2p, ■ k, (4.10) 

u = {p^- Pff = {q- kf = -2p, -Pf = -Q^ - s~i. 
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The physical kinematical region is s > M^, t < 0, u < 0. Introducing the appropriate 
kinematical factors! then, after a lengthy computation, the cross section reduces to the 
relatively simple form 



167r2 



2 



iLps^l (1 ^ ) ( 1 ^ ) + - <^£) U 



(4.11) 



where 



n=0 
oo 

- E (l^)'" [^1 (i(2^" - ^^')) + i^<^2 (i(2S - M^))] , 

n=0 



(4.12) 



and similarly for (ff and (^£. 

The cross section in terms of scaling variables is obtained by expressing the Mandel- 
stam invariants in terms of x, Q^, and the center-of-mass scattering angle 9. In the limit 
of vanishing quark mass (but retaining of course the meson's mass) we have 

s = 

X 

{l-x)Q^ -xM'^ 

2x(l-x) ^ ^ 

The physical region corresponds thus to < x < Xmax, with Xmax given by eqn.(|3.9|). It 
is now easy to check that when x = Xmax, ■§ = M'^, i = —Q^ and w = 0, so cr'^q^ vanishes 
identically, as we claimed following eqn.( |3.8| ). Integrating a'^^^{s,t) over cos^ from —1 
to 1 then gives the cross-section cr^g^(x; InQ^) which enters the expression ( |3.8| ) for the 
anomalous dimensions, and thus ( |3.14| ) for the evolution of the Gottfried sum. 



^ The cross section is dimensionless after the extraction of a scale factor which we fix in such a 
way that in the limit in which the meson is point-like the canonically normalized sphtting function 
for emission of scalar particles is reproduced. 
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Before we can do this, however, we must decide on a particular form for the vertex 
functions. From the asymptotic behaviour ( [4.6|) , and the definitions (|4.12|) we may infer 
that 

~ ^-V— , (ps ~ ^ \' \ (4.14) 



s — J-oo S s — »-oo 



we also know that both remain finite at small s. We use these facts to justify the simple 
ansatzei ips = <^(i(2s - M^)), cp^ = ip{j{M'^ - 2i)) with 

Similar forms will be assumed for the 77 and r]' vertex functions. We also make the simpli- 
fying assumption (strictly justifiable only asymptotically) that S(p^) = f-Kf^{p^)- 

The parameters A and A give the scale of transition of the vertex functions from the 
hard (pointlike) to the soft (extended) region; they may be thought of as parameterizing 
the "size" of the bound state. The normalization condition eqn.( [4.8|) fixes the strength of 
the pion-quark coupling by fixing as a function of A. Because ma may be identified 
with the dynamical quark mass (actually the "pole" mass, since E(m|) = m^) A could 
then be fixed by requiring this to be equal to the constituent quark mass. However, given 
the large uncertainty on this quantity, this only gives us a range of acceptable values of 
A, rather than fixing it completely. There is no comparable way of fixing g.,^, so we must 
treat g.,^ and A as free parameters; the dependence of the cross-section (f4 . llj ) on these is 
however rather weak except in the large region, as we will see in detail below. 

Eqn. ( ^TTl] ) shows that the cross section is dominated by the region where either the 
i-channel or the s-channel0 intermediate particle goes close to the mass shell (which it 
can never reach because of the meson mass) — the region where either t or s is small. 
Physically, this corresponds to the emitted meson and the quark being (almost) coUinear in 
the respective diagrams fig. 3a and fig. 3b, consistent with the probabilistic interpretation. 
In this region the form factors eqn. ( [4.1^ ) are essentially pointlike if the mesons are light. 



Now, eqn. (|4.13|) shows that in the t-channel this occurs when ^ = and x is small, whereas 



^ Even though the asymptotic behaviour ( [4.6| ) may only be proven rigorously for spacelike 
we assume that it may be analytically continued to timelike < 0; we will need vertex 
functions in both spacelike and timelike regions, since although s> \M'^J<\M^ (see (|1|)). 
^ The fact that the s-channel diagrams can also be resummed and factorized just as the t- 



channel ladder diagrams is discussed for a theory of scalar gluons in ref. [ 28 1 . 
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in the s-channel it occurs when x is close to Xmax (as given by eqn. (|3.9|) ). Both conditions 
are vahd for aU Q^, even though the phase-space region (i.e., the x-range) where s and t 
are cfose to their respective minimal values gets smaller as increases. 

This implies that for moderate the cross section eqn. (|4.lH ) is not only significantly 
increasing with Q^, but also it is dominated by the collinear region where the form factors 
are close to pointlike and the dominance of ladder diagrams is valid. As increases 
the cross-section increases more slowly. Even though the bulk of the cross section still 
comes from the collinear region, its growth is effectively cut by the extended tail of the 
form factor, thus leading to anomalous dimensions which have a "higher-twist" fall-off 
as inverse powers of Q^, until eventually the nonperturbative cross-section is negligible as 
compared to the perturbative one. The cross-section ( [4.11| ) with form factors behaving 
asymptotically as (|4.14|) thus provides a smooth interpolation between the region where 
strong nonperturbative evolution takes place, and the asymptotic perturbative regionl. 

So far we have taken all the quarks to be massless; the resulting cross section ( [4.11| ) 
then diverges at small s or small t, and, as we discuss in the next section, this means that the 
evolution due to Goldstone bosons diverges in the chiral limit M — 0. This is unacceptable, 
but easily cured. Since we have effectively included radiative corrections to the two- 
loop perturbative evolution diagram of fig. 2 which turn the emitted (and unobserved) 
qq pair into a bound state (a meson) we should include also radiative corrections on the 
propagating quarki. Indeed, this is required for consistency: since we are treating the 
nonsinglet pseudoscalars as Goldstone bosons, which become massless in the chiral limit, 
chiral symmetry must have been broken dynamically, and the propagating quarks must 
have a dynamically generated self energy S(p^) (as may be seen directly from eqn.( [4.7|) ) 
of the order of several hundreds of MeV. 

The quark self energy corrections can be absorbed into the vertex functions by the 
simple replacement 

x{pf + \K k) - ^^j$^^f^x(P/ + \K k) (4.16) 



^ If instead of ( [4.15| ) we were to take the cruder forms ^{p^) = 6{Is? — p^), (p{p^) = we 
would, after a chiral rotation, recover the chiral quark model [pO[ as used in [ p!9[ . This is not an 
acceptable approximation precisely because it is the fall-off of (p{p'^), as embodied in (4.15), 
which gives the higher twist evolution at large Q^. 

^ Though not of course on the external ones. 
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in the s-channel, and similarly for the t-channel. By rewriting these new amplitudes in 
the form (|4.4|) , it is not difficult to see that the cross-section retains its simple form ( [4 .111) , 



but now with 

oo 

S -j- Zj [Sj 

- (4-17) 

n=0 

where again all of the (pf are to be evaluated at |(2s — M^), and ip£, (pi are obtained by 
replacing s with t. Although the asymptotic behaviour ([4.14| ) is unchanged, in the infrared 
we now find that, since is finite for small s, ifg falls to zero linearly with s. This is, as 
we show explicitly below, just sufficient to cure the infrared divergence by suppressing the 
emission of Goldstone bosons in the chiral limit. To incorporate this behaviour, we take 
the new ansatze 

= j^^{\{2s - M^)), = ^±^^{\{M' - 2t)), (4.18) 



with ip(jP') given by ([4.15|) , and (/?s, as before. 



On the other hand, consistency with the partonic interpretation requires that quarks 
are treated as massless or quasi-massless. Hence, the present approach will only be justified 
if the pion mass is sufficiently small that the effects of chiral symmetry breaking in the 
pseudoscalar meson sector are felt, but also sufficiently large that the inclusion of the 
quark self energy in the cross-section is not quantitatively important. This will have to be 
checked explicitly below. 



5. Computing the Evolution of the Gottfried Sum 

The dependence of the Gottfried sum is obtained by integrating the meson emission 
cross-section, eqn.( [4.1I| ), over cos^, taking the first moment, and then using the result to 
evaluate the evolution factor (|3.14| ). The full analytic expression for the integrated cross- 
section with the form factors ( |4.15|) is given in the appendix. Even though the general 
analytic form eqn.(A.l) of the integrated cross section is rather cumbersome, we may 
understand the main qualitative features of the first moment by considering various limits 
of it. 
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In the large-Q limit the first moment of the cross section for any finite M behaves as 



(5.1) 



where ko and ki are dimensionless functions of M, m^, (7^^, A and A. We thus have power 
suppression of the anomalous dimension at large Q^, i.e., "higher-twist" behaviour. This 



is of course a direct consequence of the presence of the vertex functions (pi{k,p) eqn. 
in the quark-meson coupling ( ^4.2| ), which forces a deviation from the pointlike behaviour 
of the emission vertex at large Q^. The falloff is a direct consequence of the soft 

asymptotic behaviour (f4.6| ); if we had included the logarithms in the ansatze (|4.15| ), we 
would also have found logarithmic corrections to the higher-twist behaviour. Notice that 
the anomalous dimension of the Gottfried sum vanishes more rapidly at large Q than each 
of the anomalous dimensions for 77 and tt emission does, because the former is proportional 
to the difference of the latter two. 

Furthermore, it is easy to see that the behaviour of the first moment in the large and 
small meson mass limits at fixed is widely different. Indeed, as the meson mass M 
tends to zero (the chiral limit) the splitting function reduces to that for the emission of a 
pointlike pseudoscalar particle, namely (cf eqn. (A. 5)) 



a 



7*n 
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M2^o327r2 



^ V M2 



(5.2) 



This would lead to singular behaviour in the chiral limit due to the usual collinear diver- 
gence when the quark propagator in the diagrams of fig. 3 goes on-shell. When the meson 
mass tends to zero, however, we can no longer neglect the quark self energy, which, due 
to eqn.( [4.18| ) regulates this divergence in the infrared. We discuss the chiral limit in more 
detail below; here we merely remark that eqn.( |5.20 shows that when the emitted meson 
is very light the anomalous dimension obtained from eqn.( p.l2|) is constant, and the scale 
dependence of all moments is strong. 

In the opposite limit of very large meson mass the first moment of the cross section 
vanishes as 



7*n 
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M2 



In 



+ C2 



ci + O 



(5.3) 



where ci and C2 are dimensionless functions of Q^, nid, Qm ^ and A. Heavy mesons thus 
make only a very small contribution to the evolution, as expected. The fall-off of moments 
of the cross-section at large M is so rapid because although the cross-section itself only 
falls as 1/M2, 

Xmax ~ Q"^ / Al"^ . It is interesting to notice that the leading contribution is 
proportional to (7^, because the large M expansion of the cross-section begins at order x 
when (^TT = (see eqn.(A.6)). 

Summarizing, an analysis of various limits of the anomalous dimension computed from 
the cross section for meson emission confirms the naive expectations discussed in the end 
of section 2: due to the soft vertex functions, at large (compared to all other scales, 
including M) the anomalous dimension always reduces to higher-twist behaviour; for small 
meson mass (compared to all the other scales, and in particular Q^) the first moment 
evolves strongly, while for large mass the anomalous dimension is very small, implying 
that the nonsinglet anomalous dimensions eqn.( |3.13| ), and in particular the combination 
which governs the evolution of Sg, are rather large. 

We proceed now to a detailed computation of the evolution of the Gottfried sum. 
In order to do this, we must still fix the free parameters in the cross section eqn.( [4.1I| ). 
As we have seen in the previous section, may be determined as a function of A by 



the normalization condition eqn. (|4.8|) , or the expression ( [4.9|) , with the left hand side 
determined experimentally to be 0.38 ± 0.05 GeV^|^. The resulting dependence of 
on A is shown in fig. 4, and turns out to be rather weak over a reasonable range of A; 
moreover, the two independent determinations eqn.( |4.8| ) and ( [4.9| ) are compatible within 
the given uncertainties. In what follows we use for definiteness the value of determined 



as a function of A by eqn. (|4.8|) , which (see fig. 4) provides us with a lower bound on the 
anomalous dimension and therefore a conservative estimate of the evolution of the Gottfried 
sum. A is then treated as a free parameter; if is to be of the order of the constituent 
quark mass {ma ~ ~ \^p)^ ^ must be roughly in the range 0.4 ^ A ^ 0.8 GeV. 

We are thus left with the parameter g-^, which gives the relative importance of the 
derivative coupling to the pseudoscalar coupling of the meson to the quark, eqn.([4.2|). 
The sign of g-,^ is irrelevant, as is apparent from eqn. (^4.11J ). Solutions of truncated Bethe- 



Salpeter equations seem to indicate that and (^2 sjce rather small [p6| , but fall off on the 
same scale as 4>q] according to ( [4.12| ) and ( [4.15D this suggests (/tt < 1 and A ~ A. However if 
quark self energy corrections are included, then ^0 also contributes to (p ( ([4.17| )), suggesting 
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that (7^ ~ 10. We thus consider < (^tt ^ 1 to be a reahstic range of values, and take 
A = A. This choice is actually immaterial for small values of Q^, up to of order of 
a few GeV^, since the contribution of the terms proportional to (7^ is negligible in this 
region. This is due to the fact that terms proportional to g^^ yield a negligible contribution 
to the cross section in most of phase space, and in particular vanish in the region where 
the form factors are pointlike (and, more generally, whenever their arguments are equal). 
However they have a slower fall-off at large than the rest of the cross-section (because 
the piece of Kx (eqn.(A.4)) which is proportional to g\ grows as InQ^) and thus dominate 
it asymptotically. 

We have thus computed numerically the first moment of the integrated cross-section 
eqn.(A.l) for tt'^, tt^, r\ and r]' emission, and the derivatives of these first moments with 
respect to InQ^, which give the anomalous dimension for nonsinglet evolution according 
to eqn.( |3.8|) . The results are displayed in fig. 5 (cross section) and fig. 6 (anomalous 
dimension), where the dependence on the meson mass and on the parameter A (with 
(/tt = 1) are also shown. It is apparent that the contribution to the anomalous dimension 
from TT emission is indeed significantly larger than that from 77 emission, thereby leading 
to substantial evolution of the first moment for values of between 0.1 and 100 GeV^. 
At small the cross-sections evolve very slowly; furthermore, their derivatives tend to a 
common value, thereby leading to vanishing of the anomalous dimension and flattening of 
the nonsinglet evolution. For large the cross section flattens; the precise value where 
this occurs is controlled by the parameter g-^ and may vary from ~ 10 GeV^ for (/tt = 
to g2 ^ io3 GeV2 for g^ = \. 

The relative importance of various contributions to the cross section in different 
regions, as well as the dependence on (/tt are displayed in flg. 7. It appears that t-channel 
pseudoscalar emission (flg. 3b) provides the leading contribution in the very small ~ 
0.01 — O.lGeV^ region; s-channel emission (flg. 3a) provides the dominant contribution in 
the small to intermediate ~ 0.1 — IGeV^ region; and the derivative coupling terms, 
proportional to (/tt, control the large tail. Comparing flg. 7a (tt*^ emission) with flg. 7b 
[j]' emission) demonstrates explicitly the relatively slower fall-off ( ^.3|) of the derivative 
coupling terms compared to the nonderivative ones with increasing meson mass M . 

The x-dependence of the s and t channel contributions which, upon integration, lead 
to the total cross section of flg. 7a is further displayed in flg. 8. This shows that indeed 

Although these contributions are either soft, or suppressed by ifi/md, where m is the light 
current quark mass. 
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the nonperturbative evolution is dominated by the coUinear region where the leading log- 
arithmic approximation holds and the form factors are almost pointlike. It is interesting 
to observe that in the s-channel the coUinear region is dominant despite the fact that this 
is the large-a; region, and the cross-section vanishes at a; = x^iax- 

The flattening of the evolution at small identifies naturally a scale at which physical 
observables should be connected smoothly to the quark-model values, which one expects 
to hold in the infrared regime. It is pleasing to notice that the flattening occurs at ~ 
(200 MeV)^, the natural "QCD scale" at which conflnement and chiral symmetry breaking 
are expected to occur. Moreover the consistency condition resulting from ( |3.11j ) can indeed 
be satisfled in this region (and in no other); when 7^"^ ~ 0, tti ^ 0. At this scale there is no 
dynamically generated meson component, hence the sea is expected to be symmetric and 
Sg take its quark model, valence value Sq = |- This result cannot be compared directly 
with those obtained in effective models (such as Skyrme or chiral quark soliton models) 
because in such models there is no scale dependence. 

We can thus proceed to compute the scale dependence of the Gottfried sum assuming 
that it takes the quark model value Sq = | at a reference scale Qq = 200 MeV, where the 
evolution flattens, and then using eqn.( |2.7| ). We take A = A and we let this parameter vary 
in the range discussed above, while we flx from fig. 4 (solid line), i.e. from eqn.( [4.8|) ; 
we then vary < < 1. Good agreement with the experimental data is obtained for 
A ^ 550 MeV, as shown by fig. 9, where both the dependence on A and (7^ are shown. 

Fig. 9 also displays the fiattening of the evolution at both small and large Q^. Because 
of the multiplicative character of the evolution the precise small-Q^ value Qo at which the 
condition Sg{Qo) = ^ is enforced is immaterial. Indeed, we have checked explicitly that 
if instead of taking Qo = 200 MeV, we fix Qo by requiring that the anomalous dimension 
be smaller than a threshold value (so that Qo may depend on the specific values of the 
parameters) the results are essentially unchanged. The value at which the evolution fiattens 
at large is controlled by g-n-, and is generally smaller than the value at which the cross 
section for the emission of each particle separately fiattens, as we discussed above. The 
fiattening occurs thus typically around ~ 3 GeV^ in the extreme case of (/tt = 0, and 
g2 ^ 100 GeV2 ifg^ = l. 

In any case, we predict that there is a significant scale dependence of So in the region 
of ~ 1 GeV^, which could be experimentally detectable. Furthermore, unless (/tt ^ 1, it 
appears that the asymptotic value of Sq is not yet attained in the region where present-day 
data are taken; rather, the asymptotic value can be as low as S'g(oo) ^ 0.20 if (7,^ ~ 1. 
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We can now address the issue of the self-consistency of the present approach. Firstly, 
the very fact of finding a positive value of the anomalous dimensions is of itself nontrivial; 
there is no fundamental reason why the cross section for meson emission should always 
increase with . This provides an a posteriori consistency check of the probabilistic in- 
terpretation of the splitting functions which we have assumed throughout. Furthermore 
the flattening of the cross-sections at small is another nontrivial consistency check; 
without it we could not assume that the meson component of the sea is generated entirely 
by evolution, due to the condition The flattening at large is consistent with 

the asymptotic recovery of perturbative behaviour. A study of the x and dependence 
of various contributions to the cross-section (figs 7 and 8) shows that indeed significant 
dependence is generated nonperturbatively only when the i-channel (at small Q^) and the 
s-channel evolution (at intermediate Q^) are dominated by the respective quasi-coUinear 
singularities; this is consistent with the probabilistic interpretation of the evolution equa- 
tions. 

Finally, let us consider the chiral limit of our approach. In this limit the cross section 
eqn. ( |4.11| ) diverges in the infrared because the coUinear singularities are no longer regulated 
by the meson mass. This would lead to the behaviour eqn. (|5.2| ) of the splitting function, 
and thus to an anomalous dimension which does not display higher-twist behaviour at 
large Q^, but rather coincides with that obtained in a theory with pseudoscalar pointlike 
gluons. This is untenable both theoretically and phenomenologically, and would seem 
to suggest that the chiral limit reveals an inconsistency in our approach. However, as 
discussed at the end of section 4, in order to examine the chiral limit consistently it is 
necessary to compute the meson emission cross sections and the evolution of the Gottfried 



sum by using the form eqn. (|4.18D for the vertex functions, which corresponds to including 
the self-energy corrections in the internal quark propagator. The external (parton) quarks 
are kept massless, so the Mandelstam invariants retain their form eqn. (|4.13|) . 

In order to simplify the computation, we have approximated the function ipg 

3 

eqn. (|4. 181 ) with the linearly rising function = 'pf^'^ fo^ small s, and the asymptotic 
form eqn. (|4.15|) for large s, joined by demanding continuity. An analogous approximation 
is taken for (/;£. The dependence of the cross-section on the pseudoscalar meson mass with 
and without the inclusion of the self energy correction are displayed in fig. 10. It may be 
seen that the self energy correction does indeed succeed in taming the infrared divergence, 
thus giving sensible results in the chiral limit. 
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However, the physical value of M^^ is sufficiently large (and thus in practice far from 
the chiral limit) that the self energy correction makes little difference to the actual evolution 
of the Gottfried sum. This is shown in fig. 11, where the evolution of the Gottfried sum 
with and without self-energy corrections are compared (with difi'erent values of gf^). It is 
thus apparent that the partonic approach is indeed consistent with the smoothness of the 
chiral limit, since our results for the non-perturbative evolution of the Gottfried sum are 
essentially insensitive to the dynamical quark mass generation. It is interesting to observe 
that the normalization of the cross section does depend on the propagating quark mass, 
but its derivative does not. 

We conclude that our approach provides a fully self-consistent, rather stable, and 
almost parameter-free determination of the anomalous evolution of the Gottfried sum, 
which allows us to understand the current experimental data in a natural way, and suggests 
a further Q'^ dependence of the Gottfried sum which may be observable as and when new 
data at different values of Q'^ become available. 

6. Conclusion 

In this paper we have presented an extension of the Altarelli-Parisi evolution equation 
to include bound state emission, and computed the relevant splitting functions in the 
particular case of pseudoscalar emission, which gives the only significant contributions 
in the nonsinglet channel, due to the large breaking of U{Nf) fiavor symmetry in the 
pseudoscalar meson spectrum, and the relatively low mass of the pseudo-Goldstone bosons. 
Our computational method is based on the probabilistic interpretation of the Altarelli- 
Parisi equations, and leads to a determination of the non-perturbative scale dependence of 
nonsinglet quark distributions due to the direct coupling of the quarks to the pseudoscalar 
bound states. We have then determined explicitly the relevant splitting function, we have 
computed its first moment, and we have determined the scale dependence of the first 
moment of the nonsinglet quark distribution, and hence that of the Gottfried sum. 

The overall qualitative features of our results for the splitting functions and for the 
evolution of the first moment confirm a posteriori the consistency of our approach: the 
anomalous dimensions are positive definite (consistently with the probabilistic interpreta- 
tion), and flatten both at large Q^, where they reduce to higher-twist behaviour, and at 
small Q^, around the confinement scale. The evolution of the Gottfried sum is computed 
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in terms of a single free parameter, which is however fixed independently to within a fac- 
tor of two, and turns out to be exactly such as needed to explain the recent experimental 
results, in that it joins smoothly to the quark model value at the confinement scale, and 
reproduces the observed, significantly smaller value at = 4 GeV^. The evolution has a 
consistent chiral limit, and is stable upon variation of the parameters. 

We believe that the results presented here provide an appealing and consistent reso- 
lution of the puzzle posed by the recent data on the Gottfried sum; the value of the sum 
which is measured experimentally is significantly smaller than the quark model expectation 
because non-perturbative effects lead to significant evolution in the intermediate re- 
gion. This evolution produces a flavor asymmetry in the quark sea which screens the value 
of the valence asymmetry, hence the value of the Gottfried sum. The ultimate dynamical 
reason for this screening is the breaking of the U{Nf) flavor symmetry to SU{Nf) in the 
pseudoscalar meson spectrum due to the chiral anomaly. In this sense, the evolution of 
the Gottfried sum is anomalous, just as the evolution of the flrst moment of the polar- 
ized structure function [|2l| is expected to be [|2^. Our results also provide a deflnite, 
experimentally testable prediction for the scale dependence of the Gottfried sum, and in 
particular suggest that the asymptotic value of the Gottfried sum may be signiflcantly 
smaller than the present experimental determination. 

Beyond the immediate application to the Gottfried sum, our results provide a rather 
more general technique for the study of the non-perturbative aspects of deep-inelastic 
scattering. In particular, it will be interesting to compute the evolution of the full nonsin- 
glet quark densities (not just their flrst moments) which, at scales below a few GeV, should 
be rather different from what expected on the basis of perturbative QCD alone. Indeed, 
there are now experimental indications that the ratio of structure functions / F2 evolves 
differently (i.e. signiflcantly more strongly) than expected perturbatively . Also, the 



present approach could shed new light on the old idea[^ that the sea distributions are 



generated radiatively from simple valence distributions at low scales where a quark model 
valence picture of the nucleon is presumably valid. Whereas all previous attempts to this 
were hampered by the obvious limitations of perturbative evolution, the present approach 
could provide the required nonperturbative information. In particular, it is encouraging 
to note that in contrast to the perturbative evolution the flatness of the non-perturbative 
evolution at low scales will greatly reduce the sensitivity of the flnal distributions to the 
(a priori unknown) starting scale. 
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The possibility of determining in this way the small-x behaviour of the nonsinglet 
distributions seems especially intriguing. In the singlet channel, however, the bulk of the 
evolution will be controlled by the usual perturbative anomalous dimensions, since the 
nonperturbative anomalous evolution we determine, albeit unusually large for the nonsin- 
glet channel, is always small compared to the singlet perturbative evolution. Nonetheless 
the effects which we discuss here may still be important in peculiar kinematical ranges, 
such as very small x at intermediate Q^. 
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Appendix. 



The cross-section '^(x;(5^), obtained by integrating the cross-section eqn. ([4.11| ) 
over cos 6 , is given by 



1 ml{ml + A2)2 f 2Q2(i _ ^) (4Q2 _ (3^2 - AK^)x) 



71^ 



J 71 



+ 



(M2 + 4A2)2 ((52 - M2x) [2g2 _ (M2 - 4A2)x] 
x(M2x-Q2(i_a.)) 



+ 



X 



g2(i _ x)2[2Q2(i _ a;) _ (M2 - 4A2)x]2 
1 



g2 (M2 + 4A2)3 [Q2(i _ _ M'^x] [2Q2(i _ a;) _ (m2 - 4A2) x] 



4Cl(x,Q^M^A2)ln 



(1 -x) (Q2 _m2x) 



+ 



-C2(x,Q2,M2,A2)ln 
gl (m2 + A2)2 



M2x2 

x(M2(1 + x)+4A2(1-x)) 



m. 



{ml + A'^y 



(l-x)(2Q2_ (M2-4A2)x) 
x(3M2 - 8A2) - 2(M2 - 4A2)x2 - (1 - x)(5 - 4x)Q^ 



(1 - x)[2Q2(i _ a;) _ (M2 - 4A2)x]2 

1 — X 



+ 



[M2(l + x) +4A2(l-x)] 
x(l - x)[3(M2 - 4A2)x - 2(3 - x)Q2] 
2[Q2(i -x)- M2x] [2(52(1 - x) - (M2 - 4A2)x 

, / x(M2(1 + x)+4A2(1-x)) 
X In ' 



(1 -x)(2Q2 _ (M2 -4A2)x) 
in terms of the parameters defined in section 2 ( ^.131 ), ( ^.151 ) and with 



(A.l) 



Ci(x, Q^ M^, A^) = -2M2(m2 + 4A'^fx'^[M^x -Q^{1- x)] 

+ [m2 + 4A2 + (3M2 - 4A2)x] (1 - x)Q^ [{M^ - 4A^)x - Q'(l - x)] 
C2(x, Q^, M2, A2) = 2M2(m2 - 4A2)(m2 + 4A2)2x2 [m2x - Q^(1 - x)] 

- [M^ + 4A2 + (3M2 - 4A2)x] [M^ - 4A2 + 2Q^{1 - xf] . 

In the hmit of large eqn. (|A.1| ) gives 



(A.2) 



1 m2(m2+A2)^ 



7r^ 



J 71 



+ 



(A.3) 
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where 



4(1 -x) 2[M2(1 + 3x) +4A2(l-x)] , f M^{1 + x) + AA'^il - x) 



(M2 + 4A2)2 (M2 + 4A2)3 

<7Mmi + A2)2 (1-x) 



2M2x 



m 



{m\ + A2)2 [M2(i + x) + 4A2(1 - x)] 



4M2 ^ 2M2x[M2(1 + 3x) + 4A2(1 - x)] ( M'^{\ + x) + 4A2(1 - x) 



(M2 + 4A2)^ 



2\2 



+ 



9l + A 
(m^ + A2)2 



(l-x)(M2 + 4A2)J 
x(5 — 4x) x(3 — x) 



2M2x^ 



4(1 -x)2 2(1 -x) 



In 



X 



[M2(1 + x) +4A2(l-x)] 
2g2(l_a;) 



In the hmit of smaU meson mass M eqn. (|A.lJ ) yields 

1 m^(m^ + A^^^ 



7*n 

99 



(x;Q^ 



1 — x) In 



167r2 /2A4 ' ' 
while in the limit of large meson mass M eqn. (|A.1D reduces to 



+ 0(1) 



(A.4) 



(A.5) 



7*n 

99 



(x;Q^ 



2n2 



+ 



1 m2(m2 +A2 ) 
7r2 /2M2 

^(m|+A2)2 
(m^ + A2)2 



X 



2x 

g2(l - x)2" ^ V 1 



X 



x^ 



2 + 3x(l-x) 3(1 -xl 



(1-X2) 



In 



1 — X 
1 + X 



+ 



(A.6) 

The expressions ( |A.3| ) and (|A.6| ), however, are not very useful in practice since, for 
large Q^^ x^ax = \-M'^ jQ'^ ^0{M^ jQ^) and, for large M, x^ax = QV^' + 0(QV^^)- 
To obtain the behaviour of the first moment of (|A.1| ) thus it is not enough to take mo- 
ments of these equations; rather one must expand consistently in both l/Q^ and x, or, 
respectively, and x, which leads to rather cumbersome expressions which will not be 

given here. 



32 



References 



[1] J. Ashman et al., Phys. Lett. B206 (1988) 364; Nud. Phys. B328 (1990) 1. 
[2] P. Amaudmz et al, Phys. Rev. Lett. 66 (1991) 560. 
[3] K. Gottfried, Phys. Rev. Lett. 18 (1967) 1174 ; 

see also e.g. F. Close "An Introduction to Quarks and Partons" (Academic, London, 
1979) or E. Leader and E. Predazzi "Gauge Theories and the New Physics" (Cam- 
bridge University, Cambridge, 1982). 
[4] A. D. Martin, W. J. Stirling and R. G. Roberts, Phys. Lett. B252 (1990) 653; 

S. D. Ellis and W. J. Stirling, Phys. Lett B256 (1991) 258. 
[5] V. R. ZoUer, Phys. Lett. B279 (1992) 145; 

B. Badelek and J. Kwiecihski, Nucl. Phys. B370 (1991) 278; 
W. Melnitchouk and A. W. Thomas, Phys. Rev. D47 (1993) 3783. 
[6] G. Altarelh, Phys. Rep. 81 (1982) 1. 
[7] B. Badelek et al.. Rev. Mod. Phys. 64 (1992) 927. 
[8] See e.g. T. Sloan, G. Smadja and R. Voss, Phys. Rep. 162 (1988) 45. 
[9] S. Forte, Phys Rev. D47 (1993) 1842. 
[10] A. D. Martin, W. J. Stirling and R. G. Roberts, Phys. Rev. D47 (1993) 867; Phys. Lett. 
B306 (1993) 145; Rutherford preprints RAL-93-014 and RAL-93-027(1993); 
J. Botts et al., Phys. Lett. B304 (1993) 159. 
[11] D. A. Ross and C. T. Sachrajda, Nucl. Phys. B149 (1979) 497. 
[12] A. Signal, A. W. Schreiber and A. W. Thomas, Mod. Phys. Lett. A6 (1991) 271. 
[13] R. P. Feynman and R. D. Field, Phys. Rev. D15 (1977) 2590. 
[14] A. W. Thomas, Nucl. Phys. A532 (1991) 177. 

[15] S. Kumano, Phys. Rev. D43 (1991) 59 Phys. Rev. D43 (1991) 3067; 

S. Kumano and J. T. Londergan, Phys. Rev. D44 (1991) 717. 
[16] H. Walliser and G. Holzwarth, Phys. Lett. B302 (1993) 377. 
[17] M. Wakamatsu, Phys. Rev. D44 (1991) R2631; Phys. Rev. D46 (1992) 3762. 
[18] A. W. Schreiber et al., Phys. Rev. D45 (1992) 3069. 
[19] E. J. Eichten, I. Hinchliffe and C. Quigg, Phys. Rev. D45 (1992) 2269. 
[20] H. Georgi and A. Manohar, Nucl. Phys. B234 (1984) 189. 

[21] G. Altarelli, in "The Challenging Questions", Proc. of the 1989 Erice school, 
A. Zichichi, ed. (Plenum, New York, 1990) ; 

G. G. Ross in "Proc. of the 1989 Int. Symp. on Lepton and Photon Interactions at 
High Energy" , Stanford, 1989 (World Scientific, Singapore, 1989) ; 
G. Veneziano in "From Symmetries to Strings", proc. of the Okubofest, Rochester, 
1989, A. Das, ed. (World Scientific, Singapore, 1990) ; 
J. Ellis, Nucl. Phys. A546 (1992) 447c. 
[22] G. 't Hooft, Phys. Rep. 142 (1986) 357 and ref. therein. 

33 



[23] S. Forte, Nucl. Phys. B331 (1990) 1; 

S. Forte and E. V. Shuryak, Nucl. Phys. B357 (1991) 153. 
[24] R. D. Ball, Phys. Lett. B266 (1991) 473. 

[25] A. V. Manohar, Lectures at Lake Louise Winter Inst., Canada, 1992, San Diego 

preprint UCSD/PTH 92-10 (1992). 
[26] L. Durand and W. Putikka, Phys. Rev. D36 (1987) 2840; 

J. C. Collins and J. Qiu, Phys. Rev. D39 (1989) 1398. 
[27] S. J. Chang and P. M. Fishbane, Phys. Rev. D2 (1972) 1084; 

V. N. Gribov and L. N. Lipatov, Sov. Jour. Nucl. Phys. 15 (1972) 438, 675; 

see also Y. L. Dokshitzer, D. I. Diakonov and S. I. Troyan, Phys. Rep. 58 (1980) 269. 
[28] C. H. Llewellyn Smith, in "Facts and Prospects of Gauge Theories", proc. of the 1978 

Schladming school, P. Urban, ed. (Springer, Wien, 1978). 
[29] J. D. Sullivan, Phys. Rev. D5 (1972) 1732 ; 

A. W. Thomas, Phys. Lett. B126 (1983) 97. 
[30] R. Jackiw and K. Johnson, Phys. Rev. D8 (1973) 2386. 
[31] C. H. Llewellyn Smith Ann. Phys. (NY) 53 (1969) 521. 
[32] H. Pagels, Phys. Rev. D15 (1977) 2991; D19 (1979) 3080; D21 (1980) 2337. 
[33] H. Pagels and S. Stokar, Phys. Rev. D20 (1979) 2947. 
[34] R. D. Ball, Int. Jour. Mod. Phys.A5 (1990) 4391. 
[35] B. Holdom, Phys. Rev. D45 (1992) 2535; 

R. D. Ball and G. Ripka, in preparation. 
[36] R. D. Ball and J. Tigg, unpublished. 

[37] C. J. Bebek et al., Phys. Rev. D13 (1976) 25; D17 (1978) 1693. 
[38] P. Amaudruz et al., CERN preprint CERN-PPE/91-167 (1991). 
[39] G. Parisi and R. Petronzio, Phys. Lett. B62 (1976) 331 ; 

V.A. Novikov, M.A. Shifman, A.L Vainshtein and V.L Zakharov, Ann. Phys. 105 

(1977) 276 ; 

M. Gliick and E. Reya, Nucl. Phys. B130 (1977) 76. 



34 



Figure Captions 

1] The two-loop diagram which generates nonsinglet evolution of structure functions. 
When i = j the final state must be antisymmetrized. 

2] The diagram responsible for flavor symmetry breaking non-perturbative evolution. 

3] Deep-inelastic scattering off a quark which radiates a bound state 11: a) i-channel 
diagram; b) s-channel diagram. 

4] The constituent mass as a function of A, as given by eqn.s (|4.8|) (full line) and 
( ^?9|) (dashed line; the dotted lines correspond to the experimental uncertainty). 

5] First moment aj ^ of the cross-section eqn. (|4.lTD computed with A = A, g-,^ = 1 and 
rud fixed from the fuU line of fig. 4: a) H = 7r°; b) H = r]'; c) A = 600 MeV. 

6] Anomalous dimensions — ^'^1 ^ computed from the cross sections displayed in 
fig. 5. 

.7] Contributions of various terms in eqn. (f4.11| ) to the cross section displayed in fig. 5a 
with A = 600 MeV: full line, cross section with g^^ = 0; dot-dash line, cross section 
with g-K = dotted line, s-channel contribution (obtained setting = (/tt = in 
eqn. (|4.11|) ); dashed line, t-channel contribution (obtained setting = (7^^ = in 
eqn. (plD ). 

8] Dependence on x of the s-channel and t-channel contributions whose integral with 
respect to x is shown in fig. 7a, with several values of Q^: a) t-channel contributions 
(dotted lines of fig. 7a); b) s-channel contributions (dashed lines of fig. 7a). 

9] Scale dependence of the Gottfried sum, computed with A = A and fixed from the 
solid curve in fig. 4: a) g-j^ = l; b) solid curves, gr^ = 1; dotted curves, g^n — ^; dashed 
curves, (7^ = 0. The experimental value eqn.( p..2|) (uncorrected for shadowing) is also 
displayed. 

10] The approach to the chiral limit of the first moment of the integrated cross section: 
a) Without self-energy corrections; b) with self-energy corrections. The values of the 
parameters are the same as in fig. 5c. 

11] Detail of the scale dependence of the Gottfried sum: solid curve, as the solid curve in 
fig. 9b ((/tt = 1); dotted curve, same with dynamical mass correction; dashed curve, as 
the dashed curve in fig. 9b ((/tt = 0); dash-dot curve, same with the dynamical mass 
corrections. The values of the parameters are the same as in the curves of fig. 8 but 
with A = 500 MeV. 
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